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INTRODUCTION 


Much  progress  has  been  made  in  recent  years  in  the  development  of  a theoretical  framework 
for  the  formulation  of  variational  principles  [Sewell,  l9f>9;  Noble  and  Sowell,  1972;  Sewell,  197  3a, 
b;  Herrera,  1974;  Herrera  and  Bielak,  1976;  Sewell  and  Noble  1976].  A contribution  due  to  Sewell 
[197  3,  a]  allows  one  to  generate  dual  variational  principles  from  a single  functional  and  to  unify  the 
theory.  Herrera  [1974],  and  Herrera  and  Bielak  [1976],  developed  a formulation  of  variational  prin- 
ciples using  functional -valued  operators  which  simplifies  the  treatment  of  partial  differential 
equations  because  it  is  applicable  in  any  linear  space,  not  necessarily  normed,  or  with  an  inner 
product,  or  complete.  The  situation  has  been  summarized  by  Herrera  and  Bielak  [197  6], 

A sufficient  condition  for  an  operator  equation  to  admit  a variational  formulation  is  that  the 
operator  be  potential,  and  this  will  be  so  if  and  only  if  its  derivative  is  a symmetric  bilinear  func- 
tional. Variational  and  extremum  principles  are  interrelated  because  a necessary  condition  for  the 
existence  of  a stationary  extremum  is  the  vanishing  of  the  variation.  A sufficient  condition  for  a 
variational  principle  to  he  extremum  is  that  the  generating  functional  be  either  convex  or  concave. 

In  general,  a functional  is  neither  convex  nor  concave;  however,  if  the  linear  space  in  which  the 
functional  is  defined  can  be  decomposed  into  two  subspaces,  one  in  which  the  functional  is  convex 
and  another  in  which  it  is  concave,  the  functional  is  saddle,  so  that  Noble  and  Sev/ell's  results 
can  be  applied;  then  the  variational  principle  becomes  a pair  of  dual  principles.  Such  a decompo- 
sition can  be  expected  to  exist  under  very  general  conditions,  at  least  locally,  because  this  is  the 
case  for  surfaces  in  finite  dimensional  spaces. 

For  linear  operators  the  theory  becomes  especially  simple.  A sufficient  condition  for  a linear 
equation  to  admit  a variational  formulation  is  that  the  operator  bo  symmetric.  If  the  operator  is 
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either  positive  or  negative,  then  the  associated  functional  Is  either  convex  or  concave,  and 
consequently  the  variational  principle  is  an  extremum  one.  In  general  a symmetric  linear 
operator  is  neither  positive  nor  negative,  but  under  very  general  conditions  (e.  g.  when 


spectral  theorems  are  applicable),  the  linear  space  D in  which  the  operator  is  defined  can 
be  decomposed  into  two  subspaces,  one  in  which  the  operator  is  non-negative  and 
another  D in  which  it  Is  non-positive;  then  the  associated  functional  is  saddle  and  the 
variational  principle  generates  dual  extremum  principles  of  Sewell’s  type.  Accordingly,  it 
Is  possible  to  associate  dual  extremum  principles  to  arbitrary  equations  formulated  in  terms 
of  linear  symmetric  operators. 

The  main  advantages  of  using  functional -valued  operators  are  [Herrera  and  Blelak, 

1976]: 

I)  Problems  are  formulated  in  the  most  general  kind  of  linear  spaces,  which  are  not 
necessarily  normed,  or  with  an  Inner  product,  or  complete.  Most  work  in  this  field  has 
been  done  in  either  inner  product  spaces  [Noble  and  Sewell,  1972;  Sewell,  1969,  197  3a,  b; 

Sewell  and  Noble,  1976]  or  in  Hilbert  spaces  [Arthurs,  1970;  Robinson,  1971;  Collins,  1976; 

Brezis  and  Ekeland,  1976]  and  it  Is  generally  thought  that  this  Is  desirable,  if  not  essential, 

for  the  results  to  hold.  In  many  applications  the  introduction  of  the  Hilbert  space  structure  ! 

leads  to  unwarranted  complications.  This  is  not  required  when  functional-valued  operators 

are  used.  ^ 

II)  The  Introduction  of  superfluous  hypotheses  In  the  development  of  a theory  is  always 
Inconvenient,  because  frequently  they  needlessly  restrict  its  applicability. 

III)  The  symmetry  condition  for  the  potentialness  of  an  operator  can  be  extended  to 

1 

linear  spaces  for  which  no  inner  product  or  norm  have  to  be  defined  [Herrera,  197  4].  This  fact 
makes  it  possible  to  formulate  a theory  which  Is  rigorous  and  at  the  same  time  not  compli- 
cated. ; 

Iv)  Error  bounds  for  approximate  solutions  are  among  the  most  important  results  that  _ . 

the  theory  yields.  They  depend,  however,  on  simple  properties  which  are  Independent  of 

any  Hilbert  space  structure,  and  therefore  can  be  obtained  in  the  simple  setting  developed 
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by  Herrera.  General  results  for  such  bounds  are  given  In  Section  3 of  this  paper;  they  In- 
duce a metric  In  the  space.  By  using  these  bounds  It  Is  possible  to  define  weak  or  gener- 
alized solutions:  however,  these  will  not  be  discussed  here.  Certain  error  bounds  for  non- 
linear problems  have  been  described  by  Sewell  and  Noble  [1976]. 

The  application  of  the  theory  of  functional -valued  operators  to  partial  differential 
eguatlons  requires  one  to  express  them  In  terms  of  such  operators.  This  is  achieved  by 
repeated  use  of  Integration  by  parts  formulae.  Inthls  respect,  the  main  advantage  over 
standard  formulations  is  derived  from  the  fact  that  such  manipulations  can  be  carried  out 
more  systematically;  for  example,  the  operators  Involved  always  possess  an  adjoint.  A 
discussion  of  the  treatment  of  boundary  conditions  has  already  been  given  [Herrera,  1974; 
Herrera  and  Blelak,  1976];  and  In  another  way  by  Noble  and  Sewell  [1972  et  seq.  ] 

Due  to  its  generality,  the  usefulness  of  the  concept  of  saddle  functional  Is  great. 

In  applications  the  flexibility  of  this  concept  is  enhanced  by  the  fact  that  its  definition 
[Sewell  1969,  §2  (vli)]  depends  essentially  on  the  system  of  coordinates  used.  This  is  a 
fundamental  difference  with  respect  to  the  same  concept  when  applied  to  surfaces;  the 
saddle  property  of  a surface  Is  Independent  of  the  system  of  coordinates  used. 

In  some  applications  it  may  be  difficult  to  find  a decomposition  of  the  space  with  re- 
spect to  which  the  functional  Is  saddle.  For  linear  problems  [Herrera,  1974],  this  amounts 
to  constructing  a decomposition  (D^,  D ) of  the  space,  such  that  the  operator  is  non- 
negative  on  and  non-positive  on  D . However,  the  restrictions  imposed  by  these 
conditions  are  not  too  severe,  because  such  ^ decomposition  is  not  unique,  and  in  this  paper 
a wide  class  of  variational  principles  is  given  for  which  it  is  easily  obtained. 

For  initial  value  problems  only  stationary  and  not  extremum  principles  had  been  ob- 
tained, until  recently.  However  Herrera  [1974]  obtained  dual  extremum  principles  for  In- 
itial value  problems  that  were  later  applied  to  a large  class  of  problems  by  Herrera  and 
Blelak  [1976].  Later  Collins  [1976]  extended  these  results  to  a general  kind  of  dissipative 
system  by  considering  simultaneously  the  adjoint  equation.  In  France,  extremum  principles 
for  initial  value  problems  corresponding  to  certain  parabolic  equations  have  just  been 


obtained  [Brezls  and  Ekeland,  197(>];  they  Impose  more  severe  restrictions  on  the  boundary 
conditions  than  those  mentioned  above. 

The  procedure  used  by  Collins  for  dissipative  systems  (essentially,  the  mirror  method) 
to  generalize  Herrera's  principles  offers  Interesting  possibilities,  because  this  approach 
Is  not  limited  to  dissipative  systems.  In  this  paper  a general  procedure  for  deriving  dual 
extremum  principles  by  the  mirror  method  Is  developed;  It  Is  applicable  to  any  linear  prob- 
lem formulated  In  terms  of  a non-negative  (or  alternatively  non-positive)  functional -valued 
operator.  To  Illustrate  the  method.  It  Is  applied  to  the  heat  and  the  wave  equations.  To 
our  knowledge,  these  are  the  first  extremum  principles  for  the  initial  value  problem  associ- 
ated with  the  wave  equation,even  though  further  work  will  be  required  to  apply  them.  The 
principles  derived  for  the  heat  equation  are  different  from  those  available  up  to  now  and 
impose  less  restrictive  conditions  on  the  test  functions. 

A brief  description  of  the  theory  of  functional-valued  operators  Is  given  in  Section  2 
and  saddle  operators  are  Introduced  there.  Variational  and  extremum  principles  for  arbi- 
trary saddle  operators  are  given  In  Section  3 and  In  the  next  Section  these  results  are  used 
to  develop  a general  theory  of  non-negative  unsymmetric  operators.  Section  5 is  devoted 
to  extending  this  theory  to  problems  defined  in  affine  subspaces.  Some  applications  of  the 
theory  to  partial  differential  equations  are  given  in  Section  6.  Only  the  heat  and  wave 
equations  are  considered,  but  the  application  of  the  theory  to  other  equations  can  be  made 


In  a similar  manner. 
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2.  SADDLE  OPERATORS 

The  notion  of  saddle  functional  [Sowell,  19(>9;  Noble  and  Sewell,  1972]  has  been 
used  extensively  In  the  formulation  of  dual  extremum  principles  and  Herrera  [1974]  has  used 
It  to  develop  an  approach  adequate  for  the  treatment  of  partial  differential  equations 
[Herrera  and  Blelak,  1974,  1976],  In  this  section  a version  of  this  method  which  Is  suitable 
for  linear  problems  will  be  described. 

Let  D be  any  linear  space,  not  necessarily  metric,  or  normed,  or  with  an  inner  pro- 
duct, or  complete.  The  value  of  any  real  linear  functional  f;  D -♦  at  an  element  v t D 
Is  a real  number  to  be  denoted  by  <f,  v>  . It  Is  emphasized  that  this  Is  not  an  Inner 
product,  even  though  the  same  notation  is  often  used  for  such  [In  particular  by  Noble  and 
Sewell,  1972].  The  linear  space  of  all  such  functionals  will  be  denoted  by  D . 

Consider  a linear  operator  P:  D -*  D which  is  therefore  a functional -valued  linear 
operator.  Thus,  for  each  u € D there  Is  an  element  Pu  € D , and  the  value  of  this  func- 
tional at  any  v e D Is  (Pu,  v>  . Given  any  such  operator,  the  adjoint  operator 
* * 

P ; D -*  D , defined  by  the  condition  that 

( P*u,  v>  = ( Pv,  u>  (2.  1) 

is  satisfied  Vu  e D and  Vv  € D,  always  exists.  The  operator  P will  be  said  to  be  sym- 
metric,  or  self-adjoint,  whenever  P = P . 

Two  subspaces  Dj  C D and  D^  C D are  called  a decomposition  of  D If  Dj  and 
D_  are  linearly  independent  and  D = D,  + D . In  this  case,  given  any  u e D there  Is 

Z lb 

a unique  couple  (Uj,  u^)  belonging  to  Dj  and  D^  respectively,  such  that  u = Uj  + u^  : 

the  elements  Uj  and  are  called  proiectlons  of  u on  Dj  and  . Projections  of 

♦ 

linear  functionals  and  functional-valued  operators  are  defined  similarly.  Given  f t D 
and  P:  D -*  D , the  projections  fj  c D and  f^  c D are  such  that 

<fj,u>  = (f,Uj)  < f2,u)  = (f.u^)  (2.2a) 

»■  * 

for  every  u t D,  while  the  projections  Pj;  D -*  D and  P^;  D D satisfy  the  condi- 
tions 
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(2.  2b) 


PjU  = (Pu)j,  P^u  = (Pu)^ 


From  these  definitions  it  follows  that 


f = f,  . f2 


(2.  3a) 


P-P1+P2. 


(2.  3b) 


Definition  2.  1.  E C D be  a subspace  of  D . Then  an  operator,  P:  D - D is 

said  to  be  non-negative  on  E if 

< Pu,  u)  > 0 for  every  u e E ; (2.4) 

and  positive  on  E if  the  equality  sign  in  (2.4)  holds  only  when  u = 0 . When  E = D the 
specification  "on  E"  will  be  omitted. 

If  P is  non-negative  on  E the  non-negative  square  root  of  < Pu,u)  will  be  de- 
noted by  II  u lip  whenever  u e E . The  set 

Np=  {U€  e|||u|Ip=  0}  (2.5) 

may  contain  non-zero  vectors.  However,  Np  is  always  a linear  subspace  of  E C D . 
Lemma  2,  2.  Let  P be  non -negative  on  E . Then  Np  is  a linear  subspacc  of 


E C D . 


Proof.  It  is  required  to  prove: 


1)  u £ Np  =>  ou  c Np,  V or  £ R 
li)  u,  V £ Np  =>  u + V£  Np  . 


Property  1)  is  obvious  because 


< P(au),  ou>  = a ( Pu,  u)  . 


Property  11)  follows  from 


J < P(u4v),  u + v>  + J < P(u-v),  u - v>  = ( Pu,  u>  + < Pv,  v>  (2.  6) 

which  can  be  obtained  expanding  the  left-hand  member  of  this  equation.  Recalling  that 
each  of  the  terms  Involved  in  this  equation  is  non-negative,  the  property  is  clear. 


Definition  2.  3.  Let  P : D -*  D be  non-negative  on  E . Then  u,  v t D <ite  siiiii  to 
be  P-equlvalent  if  and  only  If  u-v  t Np  . In  this  case  one  writes; 

u«pV  . (2.7) 

In  view  of  Lemma  2.  2,  this  Is  an  equivalence  relation. 

Non-posltlve  and  negative  operators  are  defined  by  reversing  the  sense  of  the  in- 
equality (2.4).  For  this  kind  of  operator  there  are  definitions  and  results  corresponding  to 
those  already  given  for  non-negative  and  positive  operators.  It  Is  worth  noticing  that  the 

relation  a p becomes  equality  when  P Is  either  positive  or  negative  on  E . 

* 

Definition  2.  4.  Let  be  a decomposition  of  D and  P;  D -►  D be  a self- 

adjoint  operator.  Then,  P is  said  to  be  saddle  with  respect  to  Dj,  1^  P is  non- 
negatlve  on  Dj  and  non-positive  on  . It  Is  strictly  saddle  if  P is  positive  on 
and  negative  on  • 

It  must  be  observed  that  in  this  definition  the  sub-spaces  can  not  be  interchanged. 
However,  this  manner  of  introducing  them  simplifies  many  of  the  propositions  to  be  pre- 
sented. For  saddle  operators  the  bilinear  functional  (Pu^,  v^>  - (Pu^.v^)  is  symmetiic 
and  non-negative.  The  non-negative  square  root  of  ( Pu^,  u^)  - ( 

noted  by  llullpj  when  P is  strictly  saddle  l|u||p  is  positive  and  therefore  llu||p  is  a 
norm.  In  this  case  the  bilinear  functional  ( Pu  , v ) - ( Pu^,v  ) is  an  inner  product. 

i i £•  Lt 

Definition  2.  5.  Let  P;  D -*  D be  saddle  with  respect  to  the  decom positin.i  , 

D^  of  D . Then,  u,  v « D are  said  to  be  P-equivalcnt  if  and  only  if  ||u-v  ||  p = 0;  jri^ 
this  case  one  writes 


For  saddle  operators  the  set 


U«pV  . 


Np=  {u  € D 


(2.  8) 


is  a subspace  of  D . Thus  relation  (2.  8)  is  an  equivalence  relation. 

The  following  lemma  will  be  used  in  the  construction  of  dual  extremum  principles. 


Lemma  2.  6.  Let  P;  D D be  saddle  with  respect  to  the  decomposition  Dj, 
Define  the  subspaces  of  D 


D = {u  t d|  P u = 0} 
d 1 

= {u  e d|  P^u=  0}  . 


(2.9a) 
(2.  9b) 


Then  the  operator  P is  non-positive  in  D and  non-negative  in  D.  . Even  more,  if 

a b 

u c D 


( Pu,  U>  = - Hu  I 


(2.  lOa) 


If  u € D. 


< Pu,  U>  = Hu  I 


(2. 10b) 


Proof.  Assume  u c D_,  and  let  u,  and  u.,  be  its  projections, 
a 1 2 


Then 


< Pu,  Uj)  = 0 


and  therefore 


(Pu2,Ui>  = <Pu,u,>  - <Pu,,u,)  = -<Pu,,u,>  . 


1’  “1' 


1’  1' 


Hence 


<Pu,u>  = <Pu,u,)  + <PUj,U^>  + <Pu,,u,) 


2’  2' 


^^2’^2^  -<PUi.Ul>  = -llu||p<0  . 


The  proof  of  the  other  part  of  the  lemma  is  similar. 

In  some  applications  it  is  relevant  to  consider  operators  which  are  saddle  on  a sub- 
space  only.  Let  E C D be  a subspace,  and  Ej,  E^  a decomposition  of  E;  then  P;  D -♦  D 
is  said  to  be  saddle  on  E with  respect  to  the  decomposition  E^,  E^,  if  P is  non-negative 
on  Ej  and  non -positive  on  E^  . 

Strictly  saddle  on  E is  defined  correspondingly. 
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5.  VARIATIONAL  AND  LXTRLMUM  PRINCIPLrS 


'I 


In  this  section  a linear  space  6 will  be  considered  whose  elements  will  be  repre- 
sented by  u . Following  Herrera  [1974],  the  concept  of  derivative  of  a functional  will  be 
used  In  the  sense  of  additive  Gateaux  variation. 

A 1 /V  A 

Definition  3.  1.  Let  X:  D -*  R be  an  arbitrary  functional  and  u e D . Then 
X'(u)  t D*  is  said  to  be  the  derivative  of  X at  u If 

^ X(G  + \0)|  = <X'(u),0)  (3.  1) 

a X X=U 

/V  A 

for  every  v € D . 

The  partial  derivatives  can  be  defined  In  terms  of  the  total  derivative. 

Definition  3.  2.  Let  Dj,  D^^  be  a decomposition  of  D . Assume  the  derivative  of 
the  functional  X;  D - R^  exists  at  u € D . Then  the  partial  derivatives  X,  j(u)  € D and_ 
X,  ji(u)t  D are  defined  as  the  projection  of  X'(u)  c D , 1.  e.  in  terms  of  (2.  2) 

X,  j(G)  = {X'(G))j  and  X,  jj(G)  = {X'(G)}jj  • (3-2) 

Associated  with  a linear  functional -valued  operator  P:  D - D and  for  a given 

A A ♦ 

f t D , consider  the  equation 

pG  = f . (3.  3) 

When  P Is  self-adjoint  it  is  possible  to  establish  a variational  principle.  To  this  end, 
define  the  functional  X ; D -►  R’  by  • 

. 1 A A A A A 

X(u)  = j ( Pu,  u>  - < f,  u)  . (3.  4) 

When  X is  given  by  { 3.  4),  the  derivative  X'(u)  always  exists  and  is 

X'(S)  = pG  - f . (3.  5) 

In  view  of  this  equation,  the  following  theorem  is  obvious. 

Theorem  3.  3.  Let  P:  D — D*  be  self-adjoint.  Then  u t D satisfies  (3.  3)  if  and 
only  if 

X'(u)  = 0 . (3.  6) 


j 


! i 
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However,  this  is  not  on  extremum  prlnrtple.  \Vlien  there  Is  avalhible  a cU  ■ ’imposi- 
tion Dj,  Djj  of  D such  that  P is  saddle  vdth  n'spect  to  Pj,  Djj,  it  is  prissible  to 
foiniulatc  variational  pilnclples  that  arc  dual  and  extremum.  I'ollowing  Noble  and  Sewell 
(1972,  1976]  one  considers  the  equations  obtained  by  setting  the  partial  derivatives  of  X 
with  respect  to  and  equal  to  zero.  Tor  the  linear  equation  (3.  3),  this  leads  to 


the  system  of  equations 


A>.  A. 


X,  j(u)  = PjU  - fj  = 0 


X,  „(u)  = P„u  - f„  = 0 


(3.7a) 

(3.7b) 


These  equations,  which  are  obviously  equivalent  to  (3.  3),  can  also  be  derived  from  the 

/S  A> 

latter  equation  by  taking  projections  on  and  . 

/N  A A ^ 

Theorem  3.  -1 . Let  P;  D D be  saddle  with  respect  to  die  decomposition  D^, 

A 

of  D . Define  tho  affine  subspaces 


= {ut  d|  u satisfies  (3.7a)) 


A f A 

D . = {u  £ D I u satisfies  ( 3.  7b)  } 
b 


(3.  8a) 
(3.  8b) 


Then ; 


1)  u £ D satisfies  (3.  3)  if  and  only  if  u r p_  0 P . 

- — ■ 3 D 


li)  For  every  u ( V and  u,  c V. 

' L a a b b 


2[X(G,^)  - X(Gg)]  = llujj  - UglI-  > 0 . 


(3.  9) 


ill)  If  the  maxlir.um  of  X on  p and  the  minimum  on  p ^ coincide  and  are  attained 

— • a — ‘ ~ D " ' 

A A 

^ U,  and  U,  respectively,  then 


a P b 


(3.10) 


iv)  If_a  solution  U £ D of  (3.  3)  exists,  then 

2[X(Uj^)  - X(G^)1  = llu^  - Ull^  f llu^  - Ull-  >0  . 


(3.  11) 


•“b'  ' '“a"  " b " P " a 

In  this  r-’Tse  the  maximum  of  X on  /)  is  attained  at  u^  if  and  only  it  u^  U,  and  Hie 


minimum  of  X on  D.  is  attained  at  u.  if  and  only  if  u,  =:i.U  . In  particular  If  U and 
— b b b P 

0 are  two  solutions  of  (3.  3),  they  are  ^ -equivalent. 

Proof.  Property  i)  is  obvious. 

It  is  instructive  to  prove  part  tv)  before  part  11),  because  its  proof  is  simpler.  If 
U € D satisfies  (3.  3),  and  u^^  ( then 


A A AAA 


2X(U)  = ( PU,U>  - 2<  f,U>  = -<  PU,U> 


A AA  A AA  A , 

2X(U|^)  = ( Pu^^,Uj^>  - 2<  PU,U|^>  . 


2[X(Gj^)  - X(U)]  = < h%-V),  ' ^> 


because  P is  symmetric.  Now  u^^  = U e Dj^  (where  is  defined  by  an  extension  of 

(2.9b)),  in  view  of  the  fact  that  U as  v^ell  as  u^^  satisfy  (3.  7b):  thus,  application  of 
Lemma  2.  5 (equation  2.  lOb)  leads  to 


2[X(uJ  - X(U)J  = lluj^  - ull-  > 0 . 


Similarly 


2[X(U)  - X(G^)]  = lIGg  - U||e  > 0 . 


Adding  these  inequalities  yields  (3.  11). 

A A ^ A 

To  prove  property  ii),  let  and  then 

A AA  A A A A A A A A A 

2X(u^)  = < Pu^  - 2f,  u^>  = 2<  Pu^  - f,  u^>  - < Pu^,  u^> 

= '^all^  ■ < ^a*  = < ’^all>  ^II>  ' < ^al>  ^I^  ' ^II  ^ 


all’  air 


2 2 

It  ^ II  AA  A , II  ^ II  AA  A AA  A 

= “.n>  I’ 

where  use  has  been  made  of  the  norm  given  before  Definition  2.  5,  and  of  equations  (3.7) 


and  (3.  8).  Similarly 


2X(V  = - 2<  PC^^,  ) - 2<  u^p  . 


(3.  ]3b) 


Relation  (3.  9)  Is  obtained  by  subtracting  (3,  13a)  from  (3. 13b). 


Property  ili)  is  obviously  implied  by  li). 


4.  GENERAL  THEORY  OF  NON -NEGATIVE  UNSYMMETRIC  OPERATORS 


Let  P-  D - D*  be  a non-negative  and  not  necessarily  symmetric  operator.  Assume 

f,  g « D*  are  given  functionals,  and  consider  the  system  of  equations 

PUj  = f (4.  la) 

p*u^  = g <4.  lb) 

where  Uj.  « E)  . This  system  of  equations  may  be  written  In  terms  of  a symmetric  oper- 
ator P:  D -*  D*  where  D = D 0D  is  the  outer  sum  of  D with  Itself,  1.  e.  It  Is  the  space 
whose  elements  u = (Uj,  u^lt  D are  ordered  pairs  of  elements  Uj,  u^  c D . To  express 
equations  (4.  1)  in  the  form 

Pu  = f (4-  2) 

let  P:  D - D*  and  f € D*  be  such  that  for  every  u,  v t D 

<K1,0)  = < PUj,  v^)  + ( P^u^,  Vj)  (4.3) 

and 

<f,v>  = <f,V2>  + <g,Vj>  . (4.4) 

The  operator  P Is  clearly  self-adjoint  because 

<K1,0>  = < PUj,v^>+  < PVj,U2>  . (4.5) 

Due  to  this  fact  it  is  possible  to  formulate  variational  principles  for  the  system  (4.  1).  The 
procedure  just  described  is  a rigorous  generalization  of  the  mirror  method  [Morse  and 
Feshbach,  1953]. 

A 1 A A V- 

Theorem  4.  1.  Let  X:  D-  R be  given  for  every  u c D ^ 

X(G)  = X(Uj,U2)  = < PUp-jp  - < 

Then  equations  (4.  1)  are  satisfied  if  and  only  if 

X'(u)  = 0 . (4.7) 


Proof.  It  follows  from  Thf-orcm  3.  3,  because 

^ 1 A A A A 

X(u)  =1  ( Pu,u>  - <f,u>  . 


(4.  fi) 
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Tho  operator  P does  not  have  the  property  of  being  non-negative  or  non-posltlve 


<K1,G>  = 2(PUj,U2> 


In  general  when  a self-adjoint  operator  which  Is  functional -valued  and  linear  Is  given,  It 
Is  possible  to  find  decompositions  and  of  the  space,  such  that  the  operator  is 

saddle  with  respect  to  them  [Herrera  and  Bielak,  1976].  Finding  such  decompositions  may 
be  difficult  in  some  instances;  however,  for  the  operator  P here  considered  this  is  achieved 


easily. 


Lemma  4.2.  Let  a,  p,  y,  6 be  real  numbers  such  that  oa-Py  ^ 0,  ap  < 0,  and  -yP  > 0 . 


Define 


Dj  = {ut  D I aUj  + pu^  = 0} 


(4. 10a) 


Djj  = {ut  D 1 YUj  + 6u^  = 0} 


(4.  lOb) 


1)  The  subspaces  D^,  constitute  a decomposition  of  D . The  projections  u^ 


A,  ^ 

and  Ujj  of  any  u t D are 


Uj  = M-PvUj  - PSu^,  + aSu^) 


(4.11a) 


Ujj  = X.(q'6Uj+  P6u^,  -ovUj  - P\U2) 


(4. 11b) 


where 


or  6 -Py 


(4.12) 


11)  The  operator  P is  saddle  with  respect  to  Dj,  . Even  more,  If  P Is  positive 
then  P is  strictly  saddle. 


ill)  Let 


Y,  = -P/ff  >0,  y = y/b  > 0 


(4. 13a) 


m = — >0,  m = — 

' KV.  + -) 


(4.13b) 


. .W 


Then 


|u|l|  = mjlujllp  + m^llu^llp  . 


(4.14) 


Thus,  the  relation  u » ^ v holds,  if  and  only  If  the  relations  Uj  = p and  = p ''2 
are  satisfied  simultaneously. 

Iv)  Each  of  the  equations  (3.7a)  and  (3.7b)  is  equivalent  to 


aP  - pP  u-  = af  - pg 

Uj  2 

* 

^PUj  - 6P 


(4.15a) 


(4.  15b) 


respectively. 


Fig.l  Decomposition  of  D to  maOte  P a saddle  operator. 


Proof.  The  decomposition  of  D is  shown  schematically  in  Figure  1.  The  definitions 
(4. 10)  of  Dj  and  Djj  imply  that  these  two  sets  are  subspaces  of  D . To  prove  that  they 

A A A 

constitute  a decomposition  of  D,  it  is  only  necessary  to  show  that  u^  and  u^^  as  given 
by  (4. 11)  have  the  following  properties;  a)  t Dj  and  Ujj  i Djjih)  u = Uj  + Ujj;  c)  this 


I 

I 


Is  the  only  couple  of  elements  of  D that  has  properties  a)  and  b)  . The  proof  of  a),  b) 

and  c)  is  straightforward. 

In  view  of  (4. 9)  and  (4.  10a),  for  every  u « Dj  one  has 

^2  Z . 

-ap<Pu,u>  =a  (PUj.Uj)  +p  (Pu2.“2^  ’ 

Similarly,  for  every  u t Djj 

y6<Kl,G>  = PUj,Uj>  - Pu2’“2^  - ^ ' 

From  this  property  li)  follows  easily. 

Property  ill)  is  obtained  by  direct  computations. 

The  definition  of  projections  Implies  that  equation  (3.7a)  is  satisfied  if  and  only  if 

CPG-f,0>=O  (4.16) 

for  every  v t Dj  . Using  equations  (4.  3),  (4.4)  and  (4.10a),  it  is  seen  that  (4.16)  is 
equivalent  to 

- a(PUj  - f,  v^>  + P<P*u^  - g,  Vj>  = 0 (4. 17) 

for  every  v^  « D . Hence  (4.15a).  The  proof  of  (4.  I5b)  is  similar. 

In  view  of  this  lemma.  Theorem  3.4  can  be  applied  to  this  case. 

Theorem  4.3.  Let  X ; D -*  R be  given  by  (4.6),  a 6 - py  ^ 0,  ap  < 0 and  y6  > 0 . 
Define  the  affine  subspaces; 

D ={ue  d|  u satisfies  (4. 15a) } (4.18a) 

d 

D.  = {u  c 61  u satisfies  (4.  15b) } . (4.18b) 

b 

Then; 

1)  u = (Uj,  u^)  « D is  a solution  of  equations  (4.  1)  if  and  only  if  ‘ Pg  ^ • 

li)  For  every  = (u^,,  u^^)  r and^very  G^^  = ‘ K 

2(X(u^^)  - X(u^) ] = mj  II u^j  - Uj^j  II  ^ + m^  II  u^^  - 1 0 . (4.19) 

A 

ill)  If  the  maximum  of  X on  D_  and  the  minimum  on  /),  coincide  and  are  attained  at 

" 3 ■ O ' 

A A 

u_  and  u,  respectively 

Q O ~ 
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and  simultaneously 


al  P bl 


U u • 

a2  P b2 


Iv)  If  a solution  U = (Uj,  U^)  t D of  (4. 1)  exists,  then 

2[X(Gj^)  - X(G^))  = II  “al  ‘ ^Jlp^ 


+ ( llu^2  - ^2  II  P '*■  Il“a2  ‘ ^2  H P ^ 


(4.20b) 


(4.21) 


for  every  u t D and  u,  e D.  • 

— a a b b 

Proof.  Using  Lemma  4.2,  the  proof  of  this  theorem  is  a straightforward  application 


of  Theorem  3.4. 


5.  EXTENSION  TO  PROBLEMS  DEEINED  IN  AN  APriNE  SUBSPACE 


i 


' 


I 

i 

i 

i 

I 

i 

1, 

[ 

I 

( 


4 

X 


\ 

' \ 


In  many  applications  to  partial  differential  equations,  the  operator  Involved  may  be 
saddle  or  non-negative  when  attention  is  restricted  to  functions  satisfying  certain  boundary 
conditions.  In  such  problems  the  admissible  functions  usually  constitute  an  affine  sub- 
space and  the  theory  would  be  unduly  restricted  if  such  cases  were  not  Included. 

A set  £ C D is  said  to  be  an  affine  subspace  if  there  is  a subspace  E C D and  an 
element  w t D such  that  t = w + E . Clearly,  E = £ when  w t E,  so  that  a linear 
subspace  is  always  an  affine  subspace. 

Given  a functional  X:  D -*  r\  for  each  u t £ the  variation  of  X at  u is  a linear 

A 

functional  6X(u)  t E , such  that 

<6X(G),0>  = <X*(G),0>  (5.1) 

A A A A A A A 

for  every  v t E . When  a decomposition  Ej,  E^^  of  E is  available,  at  every  u e £ the 

A A9f: 

partial  variations  6jX(u)  € E (i  = I,  II)  are  defined  in  a manner  similar  to  partial  deriva- 
tives [Herrera  and  Bielak,  1974]. 

A 

Variational  principles  which  hold  when  the  set  of  admissible  functions  is  D Itself, 
but  which  remain  valid  when  the  set  of  admissible  functions  is  restricted  to  be  an  affine 

A A 

subspace  £ C D,  are  known  in  many  instances.  The  validity  of  such  a procedure  fre- 
quently depends  on  the  following  property  of  the  affine  subspace  [Herrera  and  Bielak, 

1976]. 

Definition  5.  1.  Let  P ; 5 -►  D*  and  f e D be  given.  An  affine  subspace  g = ^ t £ 
is  said  to  be  determinative  for  problem  (4.  2)  if  it  possesses  the  following  property; 
Whenever  u « £ , the  fact  that 

(Pu  - f,  v)  = 0 for  every  v e E (5.  2) 

implies  (4.  2). 

Lemma  5.  2.  Let  P;  D -♦  D*  be  self-adjoint,  X(u)  be  given  by  (3.  4)  and  £ C D 
be  an  affine  subspacc  determinative  for  problem  (3.  i).  Then  u « £ satisfies  (3.  3)  if  and 
only  if 
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6X(u)  = 0 . 


If  In  addition  Is  a decomposition  of  then  (5.  3)  can  be  replaced  by; 

6jX(u)  = 0 


6jjX(u)  = 0 


(5.  la) 


Proof.  Equation  (5.  3)  means 

<Kl-f,C>  = <X'(G),0>  = 0 (3.  5) 

whenever  v e E.  By  the  definition  of  determinative  affine  subspace  the  first  part  of  the 
lemma  follows.  The  other  part  is  obvious. 

The  following  result  can  now  be  easily  derived  from  Theorem  3.4. 

Theorem  5.  3.  Assume  e = w + E is  an  affine  subspace,  determinative  for  problem 
(3.  3).  Assume  in  addition  that  P:  D — D*  is  saddle  on  E with  respect  to  the  docompo- 

^ At 

sltlon  Ej,  Ejj  . The  properties  i)  to  iv)  of  Theorem  3.4  remain  valid  when  the  sets  ^ 


and  (equations  3.  8)  are  replaced  respectively  by  ; 

ft  = {G  £ ft  I u satisfies  (5.  4a) } 


(8.  6a) 


ft,  = {u  £ ft  1 u satisfies  (5.4b)} 
o 


(5.6b) 


Proof.  Under  the  assumptions  of  the  theorem  and  writing  u = w + v,  problem  (3.  3) 


is  equivalent  to 


Pv  = f - Pw 


where  v £ £,  while  P:  E -*  E*  and  f £ E*  are  the  restrictions  of  P and  f to  E . The 
theorem  follows  now  by  application  of  Theorem  3.4  to  equation  (5.  7),  after  observing  that 

|-<Pv,  v>  -<f-Pw,  v)  =X(u)  + (f,  w>  (5.8) 

That  Is,  the  functional  X(u)  differs  from  the  functional  associated  with  the  problem  (5.7) 
formulated  on  £ by  a constant  term  only. 

Theorem  5.  4.  Assume  P;  D -*  D is  non-negative  on  the  subspace  E C D . 
ft.  = w,  + E and  ft .,  = w,  + E,  with  w,,w,  t D,  be  two  affine  subspaces  of  D such 

1 1 — C c — - 1 4 ~ ■ 
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that  they  are  determinative  for  problem  s (4.1a)  and  (4 . lb)  respectively.  Define  the  affine 
subspace  e CD^  £ = Cj©£2  = ^ + ^»  where  w = ©w^  = (Wj,w^),  E = L ©E  . 

Let 

e = {u  € C I u satisfies  (5.  4a) } (5.  9a) 

d 

C.  = {u  t g I u satisfies  (3.  4b))  (5  - 9b) 

b 

Then  Theorem  4.  3 remains  true  if  the  sets  p ^ and  are  replaced  by  sets  £^  and  £ 

respectively. 

A 

Proof.  This  follows  from  Theorem  5.  3,  observing  that  the  affine  subspace  £ is  de- 
terminative for  problem  4.  2,  while  P is  saddle  on  E = E 0E  with  respect  to  the  decom- 
position Ej,  Ejj  of  E . Here 

Ej  = {0  = (Vj,  v^)  e e|  oVj  + pv^  = 0}  (5.  lOa) 

and 

gjj  = {0  = (Vj,  v^)  € eIyVj  + = 0}  . (5.  10b) 


6.  APPLICATIONS 


In  the  applications  to  be  considered,  G will  be  a region  in  n-dlmenslonal  euclidean 

space,  with  closure  G and  boundary  S = Sj  S^,  where  Sj  and  are  disjoint.  The 

linear  space  D of  functions  will  be  assumed  to  be  such  that  the  differential  operators  and 

integrals  are  well  defined  (possibly  In  a generalized  sense).  Such  would  be  the  case  if 
2 — 

the  functions  are  C on  G X[0,  T];  this  is  however,  unnecessarily  restrictive. 


A.  The  heat  equation 

For  every  u,  v e D,  the  operator  P;  D -*  D will  be  defined  by; 
T 2 

r it  9v  8u  9 u , r 9v  , 

< Pu,  V>  = / {/  ^(^  - ^ ^ + 

0 G i 1 


r 9u  9v  . 1 , /•  , 9u  9v  , , 


(6-1) 


Using  Integration  by  parts,  a convenient  expression  for  the  adjoint  operator  Is  obtained; 


9v , 9u  . 9^u  . , ^ 9u  9v 

+ ^-^)dx  - / ^ ^ dx  - 

^1 


* r ! r , 

<P  u.v)  = J^{J^ 


r 9u  9v  , 1 , r , 9u  9v  , 
h G I . 


(6.2) 


Therefore; 


f t 9u  2 , 1 r fr  1 

< Pu,  u)  = / / (^)  dx  dt  + - J {[g^  9}P\=0  ^ 


0 G 


i i 


^Wdx=llu|i; 


Assume  f , f and  f are  continuous  functions  defined  on  G X[0,  T],  S.  X[0,  T] 
G’  SI  ^2  A 


S2 


and  Sj,  X[0,T]  respectively,  and  let  f^  be  C on  G . Define  f t D , for  every  v t D, 


(6.4) 


< = /’  f 'si  S ^ / 'S2  W "2 


&v 


Bv 


0 G 


G 1 I 


Let  g , g , g and  g be  another  set  of  functions  satisfying  the  same  hypotheses  and 
G Si  S2  T 

♦ 

define  g < D by: 


(g.v)  =/■{/  "S2 


dv 


9v 


0 G 


(6.  5) 


Application  of  the  so-called  "fundamental  lemma  of  the  calculus  of  variations”  shows 


that  for  functions  € D,  the  equation 


PUj  = f 


(6.  6) 


is  equivalent  to  the  system: 


9u  9 u _ 

— - ^ ^ = r on  G X[0,  T) 
9t  9Xj  9xj  G ‘ ^ 


(6.  7a) 


9u 


(6.  7b) 


9u 

— =fs2 

0n 


(6.7c) 


!!i  !!o- 

9x,-  9xj 


on  G at  t = 0 


(6.7d) 


Similarly,  for  functions  u < D,  the  equation 


P u^  = g 


(6.  8) 


Is  equivalent  to  the  problem 
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(6.  9a) 


Tif  *’si 


■ ^G 

8t 

= ^S1 

8u^ 

8n 

" ^S2 

^"2 

9g^ 

- axj 

on  G at  t = T 


(6.  9b) 


(6.  9c) 


(6.  9ci) 


A more  standard  form  of  tfiese  problems  for  tlie  heat  equation  is  obtained  by  replacing 
equations  (6.  7d)  and  (6.  9d)  by 


Uj  = fp  on  G at  t = 0 


(6.  10) 


u=gonGatt  = T (u.ll) 

2 T 

respectively.  VWien  these  equations  hold,  the  corresponding  problems  will  be  called 
initial  and  terminal  value  problems  respectively.  Observe  that  a solution  of  problem  (6.  7) 
or  (6.  9)  may  differ  from  those  of  the  latter  by  at  most  a constant. 

The  space  D and  the  operator  P;  D - D*  is  defined  as  in  Section  4.  The  functional 

X given  by  equation  (4.  6)  is; 

T 3u  8u  a^u  9Uj 

X(Uj,u  )=/  {/  [-ar^“aT  ■ 9x.9x"  ■ y ■ ^G 

0 G 1 i 

9u,  9u  9u  9Uj  9u^ 

+ / ■ ^si’“9^r ' ^siirri'^-  ^ ^ ’ 

SI 

8u  8uj  9fp  9u^ 

- ^S2  -W  a;?;-  - 8^’-^  ’t  = 0 


8g  8u 


(6.12) 


Theorem  6,  1.  The  functions  are  solutions  of  problems  (6.  7)  ond  (6.  9) 


If  and  only  if 


X'(Uj,  u^)  = 0 . 


There  are  several  alternative  ways  of  stating  the  dual  extremum  principles  derived 
from  Theorem  4.  3.  One  that  we  have  found  convenient  is  the  following: 

Theorem  6.2.  For  every  v t D,  define; 


‘ , a^v 


“al*2’ ■ ^0*-'  ^ 2^  ^ * ax'a^  * ^ (6.13a 


0 1 1 


(6. I3t 


Let  the  set  of  functions  A C D be  such  that  v t A Implies  u , « D and  in  addition-. 


g2  T T _ 

/ v(x,  t)dt  = v{T)  + v(0)  - 2afg  + 20g.j,  + / (pg^^  - fff^)dt;  xt  G 

^ = “fsi  " P^si  ^1 


(6. 14t 


(6.  14c 


Assume  the  functional  n : A — R is  given  by; 

a 2 


» , r '^la  2 1 r "“la  "“la  . , 1 r "'o  , , 


G I i i 


G i i 


f / </  'g  ^ ‘ 4 'si  an"  ■'i  • 4 ^ 


r 0 av 

G "^Juo  ~ 


(6.  15] 


The  definitions  of  ® ^ ^ flUl  ^ . are  obtained  replacing  o-  ^ 

y and  p ^y  6 in  the  above  definitions. 


Let  vt  A and  wc  B.  Write  functions  associated  with  v by 

equations  (6.13)^^  “br“b2  functions  assoic.rted  with  w by  the  modified  ver- 

slon  of  those  equations.  Then; 


11  u . = u.  . and  simultaneously  u . s u,  if  and  only  If  u ,,u  , are  solutions 
al  bl  ' - — " ' " “ ' ' — ■ ■ ■ 


a2  bZ 


al’  aZ 


respectively  of  the  initial  and  terminal  value  problems  for  the  heat  equation. 


Z(n^(w)  -na(v)]  = mj||u^j-u^j||^p  + m2llu^2 


(6.16) 


where  ||  ||  p is  given  by  (6.  3). 


Assume  the  maximum  of  and  the  minimum  of  12^^  coincide  and  are  attained  at  v 


and  w respectively.  Then 


u , = u,  , and  u , = u,  , . 
al  bl  aZ  bZ 


(6.17) 


Thus,  they  are  solutions  of  the  initial  and  terminal  value  problems  for  the  heat  equation. 
If  solutions  Uj,  t D of  the  initial  and  terminal  value  problems  respectively  ex- 


1st,  then 


Z[l2j^(w)  - S2^(v)]  = (llubl  -^I'lp  + ll'*al  ‘^l"p^ 


.mz{l|Ubz-Uzll"p+llu^Z-^z”p>^°- 


(6.18) 


Proof.  To  apply  Theorem  4.  3,  observe  that  in  view  of  equations  (6. 1),  (6.  Z),  (6.4) 
and  (6.  5),  equation  (4.  15a)  is  equivalent  to: 

9u  9u  9^u  _ 

• -Sr  - P T - P ayl;  ■ “'g  ■ 


9Ui  9u^ 
" 9t 


“^Sl  ^ ^^Sl  SjX[0,  T] 


9u.  9u 
“-9^+  P-9^ 


!!l  !!l 

9x^  - 9Xj 


on  G at  t = 0 


9u  9g„  _ 

i = ~ on  G at  t = T 

9x.  9x, 


V.^ 


Let  V € A;  then  a straighKorwfird  computation  shows  that  u , u < D as  given  by  (6.  1 J) 

al  ac 

satisfy  equations  (6.  I9a-d).  On  the  other  hand,  equation  (6.  I9e)  is  also  satisfied  by 
u , by  virtue  of  (6. 11a). 

In  a similar  fashion,  it  is  shown  that  given  w € B,  and  satisfy  equations 

(4.15b).  Assertion  i)  follows  from  Theorem  4.  3 and  the  fact  that  the  u functions  satisfy 
equations  (6.  10)  and  (6.  11). 

When  equation  (4.15a)  is  satisfied  the  functional  X,  given  by  (4.6),  becomes 

“2’  " p ^ ^1’ “l^  ■ p ^ "'^l  ^ ■ (6.20a) 

Similarly,  when  equation  (4.15b)  Is  satisfied 

" 6 ^ ^ ^^^“2^  ■ (6.20b) 


In  view  of  these  facts  and  property  i).  Theorem  4.  3 yields  the  rest  of  the  theorem. 

Dual  variational  principles  for  diffusion  equations  were  first  obtained  by  Herrera 
[1974,  Herrera  and  Bielafc  1976].  Later  Collins  [1976]  extended  Herrera's  results  to  a more 
general  class  of  dissipative  systems.  Independently  Brezis  and  Ekcland  [1976]  have  ob- 
tained such  principles  for  a class  of  parabolic  equations.  For  the  heat  equation.  Theorem 
6.2  represents  a definite  improvement  over  previous  results  because  these  require  solving 
Laplace's  equation  at  every  time  t . Brezis  and  Ekeland's  [1976]  principles  are  in  addition 
restricted  to  the  case  where  is  void  and  u vanishes  identically  on  S . 

Dual  variational  principles  for  the  heat  equation  that  have  been  derived  in  the  past 
[Herrera,  1974;  Herrera  and  Bielak,  1976;  Collins,  1976;  Brezis  and  Ekeland,  1976]  can  also 
be  obtained  using  the  general  theory  presented  here.  For  this  purpose,  the  operator 


P:  D -»  D may  be  defined  by 
T 


.2. 


(PU,v)  -f  {/  V(^  -3^)dX4/  U^dx./  v-“dx3dt  + /v^ 

0 G 1 i SI  S2  G 


dx  (6.  21) 
t=0 

which  is  positive  on  the  linear  subspace  of  functions  with  vanishing  normal  derivatives  on 


. The  extension  of  the  theory  to  problems  formulated  In  affine  subspaces,  developed  In 
Section  5,  can  then  be  applied  to  this  operator  to  obtain  those  results.  However,  the  de- 
tails will  not  be  carried  out  here. 


B.  The  wave  eouatlon 


For  every  u,v  € D,  the  operator  P;  D -*  D will  be  defined  by: 

Bu  8v  . r 9v 


T ? 2 

^ 8v,9  u 9 u 


r , r ov  d u d u . r du  ov  . r ou  ov  , , 


/■-9u  9v  9u  9v  , , 

G <1 


(6.22) 


Integration  by  parts  shows  that 

T 2 2 

, * , /•  , r 9v  9 u 9 u . , r 9u  9v  r 9u  9v  , , 


/•  r 9u  9v  9u  9v 

* 9t  at  ^ 9x  9x  ^=T  ■ 

G 11 


(6.23) 


The  operator  P Is  non -negative  because 


/ n.  V II  l|2  1 f /r.  ou  ou  ou  - , ou  c ou  ou  , . , , , 

( Pu,  u>  = llullp  - 2 / {[(gj  ) + 3x  3x Jt=0  9t  ^ ^ 9x,  9xJt=T  ^ ^ * 


, 9u.2  9u  9u 


9u  .2  Bu  9u 


1 1 


1 1 


In  addition  to  the  functions  f and  g Introduced  In  subsection  A,  let  f^  and  g^  bo 
continuous  functions  defined  on  G . Define  f c D for  every  v t D by: 


< f,  V>  - / {/  f ^ dx  + / fgj  dx  + J r If  dx}  dt  + / [f^  ^ + 


9v 


9v 


8v 


0 G 


SI 


S2 


»v  , , 


(6.2S) 


and  g f D by 
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(6.  26) 


Then,  applying  the  fundamental  lemma  of  the  calculus  of  variations,  ft  is  seen  that  the 
equation 

Ehij  = f (6.  27) 

for  Uj  € D Is  equivalent  to  the  problem; 


1 

b\ 

a Uj 

= f^  on  G X[0,  T] 

(6.  28a) 

8t^  ■ 

ax^axi 

%. 

! 

auj 

at 

= fgj  on  SjX[0,T] 

(6.  28b) 

r 

!!i 

3n 

= fs2  ^2 

(6. 28c) 

t 

auj 

at 

= on  G and  t = 0 

(6.  28d) 

h 

Similarly,  the  equation 

aui 

afo 

= -;r-—  on  G and  t = 0 . 

axj 

(6.  28e) 

r.j 

Is  equivalent  to  the  problem; 

p%2  = g 

(6.29) 

. 

' “2 

- n on  n Yf  0 Tl 

(6.  30a) 

i 

jr 

3t^ 

ax^axj 

- 9q  A[  U,  J J 

• 

au^ 

at 

= ggj  on  Sj  X[0,  T] 

(6.  30b) 

t 

(>n 

= ^S2  ®2 

(6.  30c) 

i- 
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(6.  30d) 


on  G and  t c T 

Bu-  9g 

— - =^on  G and  t = T . (6.  30e) 

8Xj  8Xj 

As  In  the  case  of  the  heat  equation,  the  initial  and  terminal  value  problems  are  defined  re- 
placing (6.  28e)  and  (6.  30e)  by  equations  (6. 10)  and  {6. 11)  respectively;  the  observation 
made  there  about  the  relation  between  their  solutions  also  applies  here. 

The  space  D and  the  operator  P-.  D-D*  will  now  be  defined  as  in  Section  4.  The 

functional  X given  by  equation  (4.6)  is 


XfUj.u^) 


T 8u.,  8^u,  9^  u 


1 ^ 
•JJlJsr'-TJ- 


8u, 


0 G 8t  i i 

8u.  9u 


8u  — j 

/ Car  - 'si'-a;r  * <'si -5S- 

SI 


8u,  9u,  9u 

/ iiT?- Wit* 

S2 


8u  9u  9u,  9f.  9u 

/ ^^"8T  ' V ^95T  ' 9>r’  9^^tr0 

G III 

9u.  9g  8u 

^ / f^t  8t  ^ ”9x7  a>^\=T  ■ 

G f i 


(6.  31) 


Theorem  4. 1 yields  a variational  principle  for  the  problem  of  the  wave  equation  considered 
here. 

Theorem  6.  3.  The  functions  u^.  u,  t D are  solutions  of  problems  (6.  Z7)  and  (6.  29) 
simultaneously  if  and  only  if  X'(u)  =X'(Uj,  u^)  = 0 . 

Assume  a,  p,  y,  6 are  real  numbers  satisfying  the  hypotheses  of  Lemma  4.  2,  and 
define  Dj  and  Djj  by  (4. 10).  Let 
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and  equations  (4. 15b)  If  and  only  tf; 


8^w, 


8t 


8^W2 

8^ 


= Yfg  + 6g^ 


8w 

"81“  ^ ''^^Sl  * ^^Sl 


(6.  36b) 


liTr  " ’ 


all  on  the  corresponding  domains. 

Proof.  By  direct  computation,  using  the  definitions  of  P,  f and  g,  it  is  seen  that 
the  assertion  Is  true  for  functions  satisfying  equations  (6.28d,  e)  and  (6.  30d,  e)  . The 
lemma  follows  from  the  fact  that  this  set  of  equations  is  satisfied  whenever  equations 
(6.  35)  hold. 


Theorem 


6.  5.  ^ 12:  D -►  R be  defined  for  every  w = (w^,  w^)  t D ^ 


I2(w)  = ^(Wj.w^)  = X(Uj,u^) 


(6.  37) 


where  D are  given  by  (6.  33).  Assume  = (w^j,  w^^)  t D ^ ^ *'^bl’ ''bz' ‘ ^ 

satisfy  (6.  35).  In  addition,  it  will  be  assumed  that  w^  satisfies  (6.  36a),  while  Wj^ 
satisfies  (6.  36b).  Then; 

1)  Uj  and  u^  satisfy  the  initial  and  terminal  value  problems  respectively,  if  and 
only  if  w,  and  w,  satisfy  equations  (6.  36). 

2{n(Wj^)  - S2(w  g) ) = inj  / [( ) + g- g- \__T.  dx 

G 1 * 


ill) 


f,  , "l“,2-“b2'  '>l“.2-“b2>  , . 

♦ ”2  J*  — 51 — ’ — >t.o  i “ • 

G i * 


f2(Wg)  = I2(Wj^) 


(6.  38) 
(6.  39) 


if  and  only  if 


8 


(6.40a) 


r 


8u  Ou  0u  , 9Ul,i  - 

^ W 3]  _ . M on  G at  t = T 

= “3t“  ' 


and  simultaneously 


8u  8u 
aZ  bZ 


"“a2  . 


8t  ^ 8t  ’ 8x, 


on  G at  t = 0 


(6.40b) 


Iv)  If  solutions  Uj.U^*  D of  the  Initial  and  terminal  value  problems  respectively, 


exist,  then 


2{n(w^)  - fi(w^)}  = mj  j [( ) +( ) + 


8Xj  8x^ 


lt_T  "'z  f -->“•*•< 

1 G 


^^^.z-^z\z.,^(%z-^z\z 


^<%z-^z>  ^t^z-^z^  ^<V-^z>  , . 

* 8Xj  8Xj  ^ 8Xj  8Xj  n=0  ~ 


(6.41) 


Proof.  TVris  theorem  is  a straightforward  application  of  Theorem  4.  3. 

Many  extremum  principles  can  be  derived  as  corollaries  of  Theorem  6.  5.  In  many 
applications  interest  is  centered  in  the  initial  value  problem.  In  this  case  one  can  profit 
from  the  arbitrariness  of  g . Wo  do  not  intend  to  be  exhaustive  in  this  respect.  However, 
the  following  theorem  is  given  as  an  example. 


Theorem  6.  6.  Define  the  functional  tt  ; D -*  R for  every  p e D, 

T 


■IPI  ■/ f * /,'S.  ^ 4's2  If 


)dt 


+ r ff  ^ |E_]  dx  ) 

J l t 8t  8x,  8x,  ^t  = 0 ~ 

G i 1 


(6.4Z) 


Define  D C D by  the  condition  p € D if  and  only  if 


8t^ 


- K 

9xj9xj 

= 0 

on 

GX[0,T]  , 

(6.  4 3a) 

8p 

8t 

= 0 

on 

SjX[0,T)  , 

(6.4  3b) 

8n 

= 0 

on 

S2  X[0,  T]  . 

(6.4  3c) 

Assume  a solution  p. . c D of  the  initial  value  problem  for  the  wave  equation  exists.  Then: 
M 

i)  For  each  p t D 


tt(p)  > Tr(Pj^) 


(6.  44) 


li)  In  addition  the  equality  sign  holds  at  some  p€  D If  and  only  if 

- — — ■ •■II.  - ■■■'  — ■ ■ . I — — - ■ ■■  7T  ■ ■ 


(6.45) 


Proof.  Choose  the  functions  g so  that  the  right  hand  members  in  (6.  36a)  vanish, 

A A 

while  = 0 . For  functions  w c D satisfying  equations  (6.  35)  and  (6.  3oa): 

fi(w)  = < PUj,  u^)  - < g,  Uj>  - < f,  u^>  = I < PUj,  Uj>  - ^ < f,  Wj> 

because  in  this  case  equation  (4.15a)  is  satisfied. 

At  the  same  time  (6.  35b)  implies 


u,  = — w,  at  t = T 
1 or  1 


(6.47) 


Using  the  initial  conditions  (6.  35a)  gives 


1 i 

<1  ...j 


•i  "i 
T 
0 G 


(6.48) 


f 1 41',  If  ♦ 5:^  l,.„  li  ♦ / I/^Ig  If  =12  */ /si  If  '‘5 " / /S2  If  ■*2''“  > 


where  ap  stands  for  w^  . Multiply  by  2(3/ar  and  eliminate  an  irrelevant  constant  to  ob- 
tain the  functional  ir(p)  given  by  (6.  42). 
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